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ABSTRACT 

Long-duration 7-ray bursts (GRBs) are widely believed to be highly- 
collimated explosions (opening angle 9 « 1-10°). As a result of this beaming 
factor, the true energy release from a GRB is usually several orders of magni- 
tude smaller than the observed isotropic value. Measuring this opening angle, 
typically inferred from an achromatic steepening in the afterglow light curve (a 
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"jet" break), has proven exceedingly difficult in the Swift era. Here we under- 
take a study of five of the brightest (in terms of the isotropic prompt 7-ray 
energy release, E Jjiso ) GRBs in the Swift era to search for jet breaks and hence 
constrain the collimation-corrected energy release. We present mult i- wavelength 
(radio through X-ray) observations of GRBs050820A, 060418, and 080319B, and 
construct afterglow models to extract the opening angle and beaming-corrected 
energy release for all three events. Together with results from previous analyses 
of GRBs 050904 and 070125, we find evidence for an achromatic jet break in all 
five events, strongly supporting the canonical picture of GRBs as collimated ex- 
plosions. The most natural explanation for the lack of observed jet breaks from 
most Swift GRBs is therefore selection effects. However, the opening angles for 
the events in our sample are larger than would be expected if all GRBs had a 
canonical energy release of ~ 10 51 erg. The total energy release we measure for 
those "hyper-energetic" (E tot > 10 52 erg) events in our sample is large enough to 
start challenging models with a magnetar as the compact central remnant. 

Subject headings: gamma-rays: bursts - - X-rays: individual (GRB 050820A; 
GRB 050904; GRB 060418; GRB 070125; GRB080319B) 



Introduction 



Accurate calorimetry is fundamental to understanding any astrophysical phenomenon. 
In the case of long-duration 7-ray bursts (GRBs)0, three measurements are required for an 
accounting of the total relativisticE] energy release: (1) E Jtiso , the isotropic energy release 
in the prompt 7-ray emission, which is inferred from the 7-ray fluence measured by the 
detecting satellite and the associated afterglow or host redshift; (2) 9, the half-opening 
angle of the bipolar conical outflow, which is inferred from the dete ction of a characteristic 



achro matic steepening in the afterglow light curve (i.e., a "jet" break: lRhoadslll999l ; ISari et al. 



1999J); and, (3) E KE , the kinetic energy of the shock powering t he afterglow emission, which 



can be inferred either via broadband afterglow modeling (e.g., iPanaitescu fc Kumar! 12002 



1 Throughout this work, we use t he term "long-duration" GRBs to refer to those events that arise from 
the core collapse of a massive star ( Wooslev 19931: Wooslev fc Bloom 20061) . despite the fact that duration 
alone is not sufficient to distinguish fr om those GRBs associated with an older stellar population (e.g., 



Donaghv et al 



2006; 



Zhang et alJl2009h . 



2 We neglect contributions from slower moving material (i.e. 
electromagnetic emission (neutrinos, gravitational radiation, etc.). 



supernova emission) as well as non- 
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Yost et al.ll2003f), or, more accurately, from late- time radio calorimetry in the non-relativistic 



phase (e.g.. lBerger et all 12004 ; iFrail et al.l 120051 ) . 



Compilations of such measurements for the first GRB afterglows suggested that the 
collimation-corrected energy release, either from the promp t 7-rays (E^) or powering the 
afterglow (Ekf,), was tig htly clustered around ~ 10 51 erg ({Frail et al.l l200ll ; iBerger et al. 
2003al : iBloom et al.l 120031 ) . This result helped to establish the connection between GRBs 



and massive stars, as core-collapse supernovae (SNe) result in a comparable output of kinetic 



e.g. 



It further motivated efforts to utilize GRBs as standardizable cand 

Dai et alibooi iFirmani et al 



energy. 

the c osmological model of our universe (e.g., 



2007), much as has be en done for Type la SNe (IRiess et al.lll998l ; iPerlmutter et al.lll999l ; see, 



es to constrain 



2006 



Schaefer 



Filippenkol 120051 for a review 



It was soon realized, however, that the most nearby (redshift z < 0.1) GRBs were 



several orders of magn itude less energetic than the typical GRB at z > 1 (IBloom et al.ll2003 



Soderberg et al.l 12004 ). Furthermore, these underluminous events appear to be significantly 
more common (in terms of volumetric rate) than their cosmological brethren. Though the 
reason for this dichotomy remains a mystery, it suggests that perhaps long-duration GRBs 
are a more diverse population than originally envisioned. 



The launch of the Swift satellite (IGehrels et al.l 120041 ) in 2004 November heralded a po- 
tential revolution in the study of GRB energetics. With its unique combination of sensitivity 
(~ 100 GRB localizations yr _1 , an order of magnitude improvement over previous satellites) 
and precise localization capabilities (~ 3' positions arrive only seconds after the GRB trig- 
ger, with ~ 3" positions delivered minutes later), Swift promised to deliver a tremendous 
increase in the number of events suitable for detailed studies of energetics. 



Furthermore, the onboard X-ray telescope (XRT; iBurrows et al.l l2005al ) has provided 
the first detailed look at X-ray afterglow evolution. Before the launch of Swift, opening 
angles were typically inferred from the optical and occasionally radio bandpasses. X-ray 
afterglows, particularly at early times, were a relatively poorly sampled phase space. The 
additional leverage provided by the X-ray regime promised to greatly simplify the task of 
distinguishing jet breaks from other predicted spectral features in afterglow light curves due 
to the achromatic nature of this hydrodynamical transition. 

Despite these advances, measuring bolometric fluences of Swift events has proven to be 
a challengin g; task. First, the limite d bandpass (15-150 keV) of the Swift Burst Alert Tele- 
scope (BAT; iBarthelmy et al.ll2005l ) captures only a fraction of the traditional 7-ray regime. 



As evidenced by Figure CD, the uncertainties associated with Swift -E 7i i SO measurements are 
significantly larger than the pre-Swift sample, due to the difficulty in extrapolating to the 
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traditional l-10 4 keV (rest frame) bolometric bandpass. We note that the Swift measure- 
ments shown in Figure [I] incorporate a Bayesian prior on the spect ral peak energy (E p ) 
based on the E p distribution measured by the BATSE instrument (see lButler et al.ll2007l for 



details). Without this constraint the measurements would be even more uncertain. 

Second, the detailed X-ray light curves provided by the Swift XRT have revealed 
a central engine capable of injecting energy into the forward shock at late times (t ^> 
A^grb L either as short-lived X-ray flares that can contain a comparable amount of energy 
to E 1 (IBurrows et al.ll2005bl ). or as extended period s of shallow decay (so-called "plateau 



phases) inconsiste nt with standard afterglow models (IFan &: Piranll2006l ; iNousek et al.ll2006 



Zhang et al.ll2006l ). While alternative interpretations for both phenomena exist, these dis- 
coveries suggest that our simplistic adiabatic picture of afterglow evolution may need to be 
revised. 

Most importantly, surprisingly few Swift afterglows have shown the characteristic achro- 
matic steepening associated with a collimated out flow. Several groups have conducted a com 



prehensive analysis o f a large sample of X-ray (IPanaitescul 120071 ; iKocevski fc Butlerl 12008 



Racusin et al.ll2008bl ) and/or optical (ILiang et al.l 120081 ) light curves, finding that at most 
only a small fraction exhibit clear evidence for collimation. Without these collimation cor- 
rections, the true energy release horn. Swift events has remained highly uncertain (e.g., 
Kocevski fc Butlerl l200a iRacusin et al~1l2008bl ). 



Here we take a different approach. To begin with, we focus only on those Swift events 
with the largest values of E li - ]SO (Fig. [I])- In the framework of a canonical GRB energy release, 
these events should have the smallest opening angles, thereby easing to some extent the ob- 
servational bias against late-time jet breaks. Alternatively, if isotropic, these extreme events 
would place the strongest constraints on the mechanism powering these explosions. Such 
high-fluence events are also more likely to be detected by other 7-ray satellites, providing 
additional coverage in the traditional 7-ray bandpass and thereby better constraining the 
prompt 7-ray energy release. 

In addition, we only consider GRB afterglows with broadband (X-ray, optical, and radio) 
coverage extending out to late times (t > 1 month). The radio bandpass is particularly 
sensitive to wide-angle jets, as the synchrotron peak frequency typically does not reach the 
radio bandpass until days or even weeks after the burst, when the X-ray and optical bands 
may be too faint to detect a jet break. Well-sampled, broadband light curves ensure accurate 
constraints on both the opening angle and the kinetic energy powering the afterglow. 

Given these constraints, we include five events in our Swift sample: GRBs 050820A, 
050904, 060418, 070125, and 080319B. This sample is not meant to be representative of the 
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Swift population as a whole. Nor, for that matter, have we included all of the Swift events 
with large E 7 ^ iso values, as most lack th e radio and late-tim e optical coverage necessary for 
afterglow modeling (e.g., GRB 061007; ISchady et al.l 120071 ). Instead, we argue that great 
insight, in particular with regard to progenitor models, can come from studies of even a 
small number of events at the extreme. 

This work is organized as follows. In § [2] we present our observations of the afterglows 
of GRB050820A, GRB 060418, and GRB080319B. We then construct broadband afterglow 
models to extract the opening angle and afterglow energy for each one in § El To complete 
our sample, we in clude analogous res ults from previous bro adband modeling of GRB 050904 
f lFrail et al.l l2006h and GRB 070125 Jchandra et al J booah . We compare the total energy 
release from these five events with the pie-Swift sample in §0J and conclude with a discussion 
of the future of GRB energetics studies in § O 

Throughout this work, we adopt a sta ndard ACDM cosm ology with Hq = 71km s _1 
Mpc" 1 , n m = 0.27, and fi A = 1 — S7 m = 0.73 JSpereel et al.ll2007h . We defin e the flux densit y 
power-law temporal and spectral decay indices a and /3 as f v oc t~ a v~P (e.g. JSari et al.lll99a ). 
All errors quoted are la (68%) confidence intervals unless otherwise noted. 



2. Observations and Data Reduction 



2.1. GRB050820A 



GRB 050820A was remarkable in two respects. First, the Swift BAT triggered on a faint 
7-ray precursor nearly 4min before the bulk of the prompt emission, enabling contempora- 
neous 7-ray, X-ray, and optical coverage during the GRB itself. Both the X-ray and (to a 
lesser extent) the optical emi ssion abruptly brighte ned in concert with the onset of the GRB, 
suggesting a common origin (jVestrand et al.ll2006l ). 



The prompt emission was o bserved by the Kon us- Wind instrument, providing spectral 
coverage from 20keV to lMeV (jCenko et al.ll2006bl ). Extrapolating the observed spectrum 
to a rest-f rame bandpass of l-10 4 MeV, we find a fluence of (6.rt g) x 10 _5 erg cm -2 . At 
z = 2.615 (IProchaska et al.ll2007l ). the total isotropic prompt energy release in this bandpass 



was E 1)iso = (9.7i?2) x 10 53 erg. 

In addition, the X-ray and particularly optical afterglow emission from GRB050820A 
was quite bright, allowing the decay to be traced out to late ti mes. T he majority of our 
observations of GRB050820A were presented by ICenko et al.l (j2006bl ). We reported the 
detection of a likely jet break at £j = 18±2 d based on late-time Hubble Space Telescope (HST) 
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Fig. 1. — Isotropic prompt 7-ray energy release (-E^iso) of GRBs with measured redshift. All 
prompt energy releases have been transformed to the rest-frame 1 keV to 10 MeV bandpass. 
The increased sensitivity of the Swift BAT results in a population with lower values of -E 7i i SO 
and larger redshifts. It is not surprising, then, that typical Swift events should h ave large (or 



even isotropic) opening angles, making jet break measurements quite difficult (IPerna et al. 



20031 ). In this work we focus on those e v ents i n the Sw ift sample wi t h the largest values 



of R 



7,iso- 



References: 



we- Swift: 


Amat 


1 ( 


2006); Swift: 


Butler et al. 


Golenetskii et al. 


( 


2009h: 


Rau et al. 


(2009h. 



observations, lat er than nearly all previously detected jet breaks in the optical bandpass 
rtZeh et alJliooS ). 



Here we supplement this already rich data set with additional late-time X-ray and 
optical imaging. Swift XRT data were taken from the online compilation of N. ButleJE 
These detections extend the X-ray coverage out to t ~ 46 d, well past the previously claimed 
jet break time. 



3 http:/ /astro. berkeley.edu/~nat/swift; see 



Butler fc Kocevski 



(120071 ) for details. 



- 8- 



We have also obtained optical imaging of the host galaxy of GRB050820A using the 
Wide Field Channel of the Advanced Camera for Surveys (ACS) on HST (Fig. [2]). Under pro- 
gram GO-10551 (PI: Kulkarni), we obtained a total of 2238 s of exposure time in the F625W 
(Sloan r 7 ) filter, 4404 s of exposure time in the F775W (Sloan i!) filter, and 14280 s of exposure 
time in the F850LP (Sloan z 1 ) filter beginning on 2006 June 5 (UT dates are used through- 
out t his paper). We processed the data using the multidrizzle routine (IFruchter fc Hook 
2002 ) in the stsdas IRAF0 package. We used pixfrac = 0.8 and pixscale = 1.0 for the 
drizzling procedure, resul ting in a pixel s cale o f 0.05" pixel -1 . Following the recipe for point- 
soured^] photometry from Sirianni et al. ( 2005 ). we measure the following (AB) magnitudes: 
F625iy = 26.04 ± 0.13; F775W = 26.09 ± 0.11; F850LP = 25.91 ± 0.11 (i ncluding a cor- 
rection for the small amount of Galactic extinction: E(B — V) = 0.044 mag; ISchlegel et al. 



19981). These r e sults are consistent with, although slightly brighter than, the values reported 
bv lChen et all (bopgh . 



The detection of the host galaxy allows us to subtract its contribution from the afterglow 
measured at t ~ 36 d with ACS. As can be seen from Figure [2], the host contribution at this 
epoch is significant and will affect the jet break time measured in § [31 

The combined X-ray, optical, and radio light curves of GRB 050820A are shown in 
Figure [61 



2.2. GRB 060418 



GRB 060418 was detected by the Swift BAT at 03:06:08 (IFalcone et al.ll2006ah. The 
ray li ght curve shows three overlapping peaks with a total duration t 90 = 52±1 s (jCummings et al 



20061 ). GRB 060418 was also bright enough to be detected w ith Konus- Wind, 



Dut no un- 



2006). We 



certainties were provided on the derived spectral parameters (jGolenetskii et a 
theref ore use the rest -frame l-10 4 keV fluen ce derived from the Swift-BAT by 
(120071 ). At z — 1.49 (IProchaska et all 120071 ). the total isotropic prompt energy release from 



Butler et al. 



GRB 060418 was E, 



7, ISO 



x 10 53 erg 



The XRT promptly slewed to the burst location and detected a fading X-ray count erpart 
at a = 15 h 45 m 42.8 s , 5 = -03°38'26"1 (J2000.0; 5'.'8 error radius; IFalcone et alTl2006ah . Like 
many Swift X-ray afterglows, the light curve exhibits a bright flare at t ~ 128 s super-posed 



4 IRAF is distributed by the National Optical Astronomy Observatory which is operated by the Associ- 
ation for Research in Astronomy Inc., under cooperative agreement with the National Science Foundation. 



3 The host galaxy is only marginally extended in the HST images. 



Fig. 2.— HST Imaging of GRB 050820A. Left: An F850LP ACS image taken on 2005 
September 26. Center: The same field in an image obtained on 2006 June 11. Right: A 
digital subtraction of the two images, revealing the residual afterglow emission. The host- 
galaxy contribution to the afterglow flux at the time of the 2005 September 26 image is 
significant (~ 35%). 



on a power-law decay (IFalcone et al.l l2006bl ). In Figure [7] we plot the X-ray light curve 
evolution for t > 0.1 d, obtained from the online catalog of N. Butler. 



The automated Palomar 60 inch (1.5 m) telescope (P60; ICenko et al.l l2006al ) began ob- 
serving the afterglow of GRB 060418 in the Vc, Re, an d $ filters beginning 2.7 hr after the 
burst (when the source became visible at Palomar Observatory). P60 data were reduced i n 
the IRAF environment using our custom real-time reduction pipeline (ICenko et al.ll2006al ). 
Where necessary, coaddition was performed using Swarp§. Afterglow magnitudes were cal- 
culated with aperture photometry using an inclusion radius roughly matched to the full 
width at half- maximum intensity (FWHM) of the stellar point-spread function (PSF). Pho- 
tometric calibration was performed relative to the calibration files provided by A. Hender0, 
resulting in root-mean square (rms) variations of < 0.05 mag in all filters. Photometric and 
instrumental errors have been added in quadrature to obtain the results presented in Tabled 

Additional optical imaging was obtained with two large ground-based facilities to sup- 
plement the P60 light curves at late times: the Large Format Camera (LFC) mounted on 
the 200 inch (5.1m) Pa lomar Hale telescope, and the Low Resolution Imaging Spectrometer 
(LRIS; lOke et al.lll995l ) mounted on the 10 m Keck I telescope. All data were reduced in a 
manner similar to the P60 images using standard IRAF routines. 



Late-time observations of GRB 060418 were obtained with the Wide-Field Camera 



6 See http://terapix.iap.fr/soft/swarp 
7 Availablc via anonymous ftp at ftp.aavso.org. 
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(WFC) channel of the ACS on HST (GO-10551; PI: Kulkarni). The images were processed 
in an identical manner to § 12.11 There is no evidence for host-galaxy emission directly coin- 
cident with the afte rglow location. How ever, several nearby sources, which may be related 
to the host galaxy (jPollack et al.ll2009l ). may contaminate the afterglow photometry. The 
results of these observations are shown in Table [7J and Figure [3j 



The 1.3 m Peters Automated Infrared Telescope (PAIRITEL; iBloom et al.ll2006l ) began 
observing the afterglow of GRB 060418 at 5:25:34 on 2006 A pril 18. Full details of the 
PAIRITEL observations are presented by iPollack et al.l (120091 ). Here we present the full 
multi-color PAIRITEL light curve of GRB 060418, derived using ap erture photometry an d 
calibrated with respect to the Two Micron All Sky Survey (2MASS; ISkrutskie et al.l 120061 ) . 
The results of our analysis are shown in Table [71 

Finally, we began observations of the fading optical counterpart of GRB 060418 with 
the Very Large Array (VLA)@ approximately 1 d after the burst. The results of this and 
subsequent monitoring for 68 d after the burst are summarized in Tabled! For the majority 
of the observations the antennae were in the A configuration; the sole exceptions are the 
data points on 2006 June 8 (BnA) and 2006 June 25 (B). All observations were reduced with 
the Astronomical Image Processing Software (AIPS) in the standard manner. 



The resulting X-ray, optical, and radio light curves of GRB 060418 are shown in Figured 



2.3. GRB080319B 



GRB080319B, or the "naked-eye burst," has been discussed extensively in the litera- 
ture. The most remarkab le aspect of this uniq ue event was the bright optical flash (peak 
optical magnitude of 5.3; iRacusin et al.l l2Q08al ) that accompanied the prompt 7-ray emis- 
sion. Because of the temporal coincidence, the contemporaneous 7-ray and optical emission 
are believed to derive from the same physical region, with the observed 7-rays being gener- 
ated by Compton sca t tering from the same relativ istic electrons that cause the optical flash 
( IRacusin et al.ll2008al : Kumar fc Panaitescull2008l ). 



The best constraints on the prompt emission come from the Konus- Wind satellite. The 
measured (20 keV - 7 MeV) fluence does not change appreciably when converti ng to the stan 
dard bandpass. The resu lting isotropic prompt energy relea se at z = 0.937 (IRacusin et al. 
2008al : lD'Elia et al.ll2009l ) is £ 7iiso = (1.44 ± 0.03) x 10 54 erg (IRacusin et alil2008ah . 



8 The VLA is operated by the National Radio Astronomy Observatory, a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, Inc. 



Fig. 3.— HST Imaging of GRB 060418. Left: An F625W ACS image taken on 2006 
May 9. Center: The same field in an image obtained on 2006 July 11. Right: A digital 
subtraction of the two images, revealing the residual afterglow emission. There is no sign 
of any host emission coincident with the afterglow; however the photometry is complicated 
by contamination from several nearby sources (which may be related to the host galaxy; 
Pollack et al~1l2009h . 



For our modeling of GRB080319B, we dra w on the rich dat a sets of iRacusin et al. 



f l2008af ). Bloom et al.l (120091 ) . iTanvir et all (12008( 1. and ICenko et all (|2009h . To supplement 
these results, we obtained target-of-opportunity observations of the afterglow of GRB 0803 19B 
in peak-up imaging mode (i.e., with the blue 15.8 /xm filter) of the Infrared Spectrograph 
(IRS) on the Spitzer Space Telescope through a Director's Discretionary Time proposal. We 
obtained 60 dithered pointings, each consisting of two 30 s cycles, beginning at 2008 March 
21.81 (~ 2.55 d after the burst). We clearly detect the afterglow at this time with a flux 
density of 35.7 ± 3.9 /iJy (ITeplitz et al.l 120081 ). A second set of observations was obtained 
beginning at 2009 March 29.89, this time using 120 pointings. No source was detected to 
a (la) limiting magnitude of 4.0/iJy. The results of these observations are summarized in 
Table H 

We also observed the location of GRB 080319B with the Combined Array for Research in 
Millimeter-wave Astronomy (CARMA) at 95 GHz on 2008 March 20 ( meantime 11:3 ). Th e 
configuration and data reduction are identical to those described by IChandra et al.l (120081 ). 
We report a non detection at the optical afterglow position with a 2a limit of 0.50 mJy 
feock et alJl2008h . 



GRB080319B w as observed with the V LA at 4.8 and 8.5 GHz at two epochs in the first 
week after the burst (jSoderberg et al.ll2008l ). Both observations took place when the array 
was in the C configuration. We also observed the afterglow location for 33 hr between 2008 
December 20 and 2009 Jan 4 at 1.46 GHz to search for late-time emission. The afterglow 
was not detected to a 2a upper limit of < 28 /iJy. 
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Table 1. Radio Observations of GRB 060418 



UT Date a 


Time Since Burst b 


Frequency 


Flux Density 




(d) 


(GHz) 


(MJy) 


2006 Apr 19.50 


1.37 


8.46 


105 ± 45 


2006 Apr 22.43 


4.30 


8.46 


148 ± 44 


2006 Apr 25.42 


7.29 


8.46 


242 ± 59 


2006 Apr 27.25 


9.12 


8.46 


113 ± 43 


2006 May 4.15 


16.02 


8.46 


184 ± 39 


2006 May 6.20 


18.07 


8.46 


158 ± 31 


2006 May 9.17 


21.04 


8.46 


< 84 


2006 May 11.17 


23.04 


8.46 


136 ± 36 


2006 Jun 8.21 


51.08 


8.46 


< 78 


2006 Jun 25.27 


68.14 


8.46 


< 62 



a UT at midpoint of exposure. 

b Time from mid-point of exposure to Swift-BAT trigger. 
c Upper limits are reported as 2<r rms per beam area. 



Table 2. Radio/Sub-mm Observations of GRB 080319B 



UT Date a 


Time Since Burst b 


Frequency 


Flux Density c 




(d) 


(GHz) 


(A*Jy) 


2008 Mar 20.48 


1.11 


95.0 


< 500 


2008 Mar 21.56 


2.30 


4.86 


204 ± 40 


2008 Mar 21.56 


2.30 


8.46 


232 ± 42 


2008 Mar 21.81 


2.55 


1.90 x 10 4 


35.7 ±3.9 


2008 Mar 26.48 


7.22 


4.86 


167 ±61 


2008 Mar 26.48 


7.22 


8.46 


< 110 


2008 Mar 29.89 


10.63 


1.90 x 10 4 


< 8.0 


2008 Dec 25.83 


281.57 


1.46 


< 28 



a UT at midpoint of exposure. 

b Time from mid-point of exposure to Swift-BAT trigger. 
c Upper limits are reported as 2a rms per beam area. 
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The results of our millimeter and radio monitoring of GRB 080319B are shown in Fig- 
ure [8] and Table [2J along with the previously published X-ray and optical light curves. 



3. Broadband Modeling Efforts 



In the standard "fireball" formulation (e.g.. |Piranll2005l ). afterglow emission is powered 
by synchrotron radiation from relativistic electrons in the circumburst medium accelerated by 
an outgoing blast wave. The resulting spectrum is well described as a series of broken power- 
laws with three characteristic frequencies: v a , the frequency below which the radiation is self- 
absorbed; u m , the characteristic frequency of the emitting electron s; and v c , the frequen cy 
above which electrons are able to cool efficiently through radiation (IGranot fc Saril 120021 ) . 



The temporal evolution of the afterglow depends on the density profile of the circum- 
burst medium. We consider here two possibilities: a constant density circumburst medium 
[p(r) oc r°], as would be expected in an environment similar to the interstellar medium 



(ISM: ISari et al.lll998l ). and a wind-like environment [p(r) oc r 2 ], as would be the case for 
a massive star progen itor shedding its outer envelope at a constant rate before core collapse 
Jchevalier fc Aoooh . 



( Rhoads 


1999; 


Sari et al. 


1999) 



At early times, observers only notice emission from a narrow cone (opening angle 9 ~ T -1 , 
where V is the Lorentz factor of the expanding shock) due to relativistic beaming. The 
resulting evolution therefore mimics an isotropic explosion. As the shock slows, however, 
lateral spreading of the jet becomes important, and the observer eventually notices "missing" 
emission from wider angles. This hydrodynamic transition manifests itself as an achromatic 
steepening in the afterglow light curve. Measuring the time of this jet break (tj) allows us 
to infer the opening angle of the outflow (9). 

Our objective here is to translate the observed three critical frequencies, together with 
the peak flux density, F UTaaiX , and the jet break time, tj, into a physical description of the 
outflow. In particular, we shall attempt to estimate seven parameters: £"ke, the kinetic 
energy of the blast wave; n, the density of the circumburst medium; e e , the fraction of the 
total energy apportioned to electrons; e#, the fraction of the total energy apportioned to 
the magnetic field; p, the electron power-law index; Ay, the host-g alaxy extinc t ion; a nd 
9, the jet opening angle. We make use of the software described by I Yost et al.l (120031 ). a 
multi-parameter fitting program incorporating the standard afte rglow formulation , as well 
as corrections for radiative losses and inverse-Compton emission (ISari fc Esinll200ll ). 



To account for differences in instrumental configurations, we have applied a 5% cross- 
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calibration uncertainty to all data points before calculating the models. All reported uncer- 
tainties have been determined using a Markov-Chain Monte-Carlo analysis with 1000 trials 
and represent only statistical errors associated with the fit. Systematic errors associated 
with model uncertainties are likely much larger and difficult to estimate. 



3.1. GRB050820A 



The optical and X-ray light curves of GRB 050820A exhibit a dramatic r ebrightening at 



t ~ 220 s, both jumping i n concert with a strong rise in the 7-ray emission (IVestrand et al. 



20061 ; ICenko et al.l l2006bl ) . We therefore remove all X-ray and optical points at early times 
(t < 0.1 d) from our fitting routines. Likewise, the radio li ght curve exhibits a bright flare 
at t « Id that is probably due to reverse shock emission (ICenko et al.ll2006bl ). Since our 
modeling software only includes contributions from the forward shock, we include only radio 
observations at t > 5 d in our models. 

The best-fit afterglow models for GRB 050820A are plotted in Figure El an d the rel evant 
physical parameters are provided in Table [3j As discussed by ICenko et al.l (l2006bl ). the 
optical light curve exhibits a distinct break between the last ground-based optical detection 
at t ~ 7d and the HST observations at t ~ 36 d. The additional late-time X-ray data firmly 
establish the presence of this break in the X-ray bandpass as well, cementing the explanation 
as a jet break. The radio data at late times are not sufficient to distinguish between a beamed 
and isotropic outflow. 

For a constant density circumburst medium, we find two solutions with somewhat dif- 
ferent model parameters but similar overall fit quality, and both are presented in Table [3J 
Though values of the electron index p less than 2 require a somewhat artificial cutoff to keep 
the total energy carried by t he electrons finite, we slightly prefer the model with p = 1.75 
(see, e.g., iBhattacharyal l200ll for a discussion of afterglows with electron index p < 2). Here 



Table 3. GRB050820A Forward Shock Best Fit Parameters 



Medium Type a 


^KE.iso 
(10 52 erg) 


n 

(cm -3 ) 




(%) 







V 


A v (host) 
(mag) 


X? (d.o.f.) 


ISM-1 
ISM-2 


C07+8O 
41018 


n is+ ' 12 
(3.9 ± 1.0) x 10~ 4 


0.13 ±0.02 

14+ 002 
U - L —0.01 


0022+ 0011 
u.uuzz 0004 

! O + 03 
L - a -0.1 


6 6+ ' 5 
3.0 ±0.1 


1.75 ±0.02 
9 on+0 03 

z - ou -0.04 


11+ 003 
u - ii -0.02 

< 0.049 b 


1.34 (76) 
1.31 (76) 



a Radial profile of the circumburst medium, p oc r° for an ISM-like medium; p oc r 2 for a wind-like medium. 
b We are only able to calculate an upper limit for the host extinction in this case. 
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the density scale is more in line with our expectation that long-duration GRBs should in- 
habit dense regions of recent star formation, as well as previous GRB afterglow modeling 



results (e.g., iPanaitescu &: Kumarll2002l ; lYost et al.ll2003l ). Though the fraction of the total 



energy partitioned to the magnetic fie ld, £r, is small, comparable values have been inferred 
for several other previous GRBs (e.g.. IPanaitescu fc Kumar! |2002| ). 



Most importantly for our purposes, the kinetic energy of the afterglow, E KE , and the jet 
opening angle, 8, are relatively similar in the two models. Though both require extremely 
large isotropic kinetic energies, the values inferred here are only an order of magnitude 
larger than the prompt 7-ray energy release. This modest 7-ray conversion efficiency (r/ 7 = 
I \E y + Ekf] ~ 1 0% ) is consistent with theoretical predictions of the internal shock model 



flKobavashi et al.lll997l : iDaigne fc Mochkovitchlll998f ) 



The jet break time occurs somewhat earlier than originally suggested by ICenko et al. 



( 12006bl ). This is due to a combination of improved constraints from the X-ray afterglow 
coupled with a shallower post-break decay index inferred for the p = 1.75 model. 

We were unable to obtain any high-quality fits to the afterglow of GRB 050820A assum- 
ing a wind-like circumburst medium. We shall return to the lack of evidence for massive-star 
progenitors from GRB afterglow modeling in § [51 



3.2. GRB 060418 



The most striking 



t « 300 s flFalcone et al 



eature in the light curve of GRB 060418 is a bright X-ray flare at 
2006bh. Such flares have been reported in a large fraction of Swift 



XRT light curves (IFalcone et al.l 120071 ). an d are widely believ ed to be caused by late-time 
energy injection from the central engine (IZhang et al.l 120061 ). These X-ray flares can in 
some cases contribute a significant fraction of the prompt energ y release to the tota l energy 
budget, and therefore have a large effect on the post-flare decay (IFalcone et al.ll2007l ). Rapid 
variability in the X-ray light curve of GRB 060418, inconsistent with standard afterglow 
models, is seen as late as several hours after the burst. Like GRB050820A, we therefore 
only include observations at t > 0.1 d in our broadband modeling analysis. 

The resulting fits and best-fit parameters are shown in Figure [7] and Table HI Again a 
clear break is seen in the optical light curve at late times (t pa 8 d) . The X-ray afterglow 
has dropped below the XRT threshold at this time. However, the radio afterglow is still 
detected and exhibits some evidence for a steepening decay consistent with that seen in the 
optical. Unfortunately, the break time is not very well constrained, either in the optical or 
radio bands. Though not entirely conclusive, we consider this relatively strong evidence in 
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support of collimation. 



Unlike most previously modeled afterglows (c.f., GRB 050904; iFrail et al.1 120061 ). our 
results indicate that the electron cooling frequency, u c , fell below the optical bands over 
the duration of our observations. The forward-shock emission above v c is independent of 
the circumburst medium profile, leading to indistinguishable fits in the X-ray and optical 
bandpasses. While the radio behavior is divergent at early and late times, our observations 
are not sufficient to conclusively distinguish between the two models. Since the wind model 
provides a slightly better fit and does not require microphysical parameters held fixed at or 
near equipartition, we shall adopt it for the remainder of the work. 

The primary drawback of the wind-like scenario, however, is the extremely high 7- 
ray efficiency. Somehow the physical process generating the prompt emission must have 
been capable of converting ~ 99% of the outgoing blast-wave energ y to 7-rays, while most 
internal shock models predic t a maximum 7-ray efficiency of ~ 10% (IKobayashi et al.lll997l ; 
Daigne fc Mochkovitchl Il998l ) . We shall return to this issue in § HJ 



3.3. GRB080319B 



Several groups (e.g.. iRacusin et al.ll2008al ; iBloom et al.ll2009t IWozniak et al.ll2009l ) have 
presented detailed observations of the early afterglow of GRB080319B, revealing a complex 
behavior not easily understood in the context of the standard fireball model. In particular 
we note that the optical spectral index, /3o, evolves dramatically at early times, and as late 
as 0.5 d af ter the burst is to o shallow to be accommodated by our forward shock models 
(Po ^ 0.2; IBloom et al.ll2009i ). We therefore consider the evolution of the afterglow only at 
t > 0.5 d. 



Table 4. GRB 060418 Forward Shock Best Fit Parameters 



Medium Type a I E KE .is 



n I A t 



Av(host) x? (d.o.f.) 





(10 52 erg) 


(cm 3 / g cm 1 ) 




(%) 


(°) 




(mag) 




ISM 
Wind 


9+1.2 
-0.6 
12 + ' 03 

u - iz -o.oi 


n 1 H +0. 08 
0.35 ± 0.12 


0.33 c 

06+ - 01 
u - UD -0.02 


31+ 2 
15li 4 


13.9jf.6- 
22.5l«; 9 5 


2.06 ±0.01 
1 07+0.02 

1 - y ' -0.04 


0.08 ± 0.03 
< 10" 4d 


1.15 (82) 
1.10 (81) 



a Radial profile of the circumburst medium, p <x r° for an ISM-like medium; p oc r -2 for a wind-like medium. 
b For a wind-like medium, the density parameter is better known as A*, where p = 5 X 10 11 A*r -2 . 
c The best fit was achieved with e e fixed to its equipartition value. 
d We could only calculate an upper limit for the host extinction. 
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In addition, the late-time (t > 10 d) optical light curve exhibits a pronounced red bump 
not seen in the X-rays that has been a ttributed to emission from an underlying supernova 
(ITanvir et all 120081 ; iBloom et all 120091 ) . Such features have now been seen in many rela- 
tively nearby GRB optical afterglows (IZeh et all 12004 ). and are not accounted for in our 
synchrotron formulation. We therefore leave all r 1 , i', and i measurements out of our models 
at t > 10 d. 

The resulting fits are plo tted in Figure El wit h the model parameters presented in 



Table El As was pointed out by lRacusin et all (j2008al ) and most convincingly by lTanvir et al. 
(120081 ). both the X-ray and optical light curves exhibit a break at t ~ 12 d. The break is 
only clearly visible in the redder optical bands after including the contribution from the 
underlying supernova. Here the radio observations are unable to provide any constraints on 
the presence of a jet break, as the radio afterglow was comparably faint and only detected 
over the first week of observations. 

We find that two models provide a reasonable fit to the broadband data, and both re- 
quire a wind-like circumburst environment. The parameters derived for the fir st model, with 
p = 2.1 0, are broadly similar to those derived for the "wide" jet (see below) by lRacusin et al. 



(j2008al ). with one notable exception: e# differs b y two orders of magnit ude (fixed at equipar- 



tition in our model, compared with 3 x 10 in iRacusin et al.ll2008af) . Taken together, we 



find that the model parameters provided in IRacusin et all (j2008al ) provide a relatively poor 
fit to the late-time data, particularly in the optical bands. 

An alternate model, with p = 1.85, provides a marginally better fit to the data. Since 
the inferred density is slightly larger and more in line with previous GRB afterglows, we 
adopt this as our preferred model for the afterglow of GRB080319B. We note again that 
the discrepancy between the two competing models is relatively small with respect to the 
afterglow energy and opening angle. 

Table 5. GRB080319B Forward Shock Best Fit Parameters 



Medium Type a I -Ere, is 



A, 



A v (host) x ? (d.o.f.) 





(10 52 erg) 


(cm 3 / g cm 1 ) 




(%) 


(°) 




(mag) 




Wind-1 
Wind-2 


8.0t°;l 
4.9l 3 ;? 


U.UUOD_q 0003 

0.015 ±0.005 


0.11 ±0.01 

0099+ - 005 
u.uuyy_ 0001 


CO CO 
CO CO 


3.6 ±0.2 

7 0+ 07 
'■ u -o.i 


2.10 ±0.02 
1.85 ±0.02 


< 0.04 d 

07+ 04 
u - u ' -0.02 


1.13 (157) 
1.12 (157) 



a Radial profile of the circumburst medium, p oc r° for an ISM-likc medium; p oc r~ 2 for a wind-like medium. 
b For a wind-like medium, the density parameter is better known as A», where p = 5 X 10 11 J 4*r -2 . 
c The best fit was achieved with eg fixed to its equipartition value. 
d We could only calculate an upper limit for the host extinction. 
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Racusin et all (j2008al ) have incorporated the early-time data into their model of GRB 0803 19B 



by invoking a double-jetted system: the high-energy emission is focused in a narrow (9 ~ 
0.2°) jet that dominates the early afterglow (t < 0.5 d), while the material at lower Lorentz 
factor powering the late-time afterglow is channeled into a wider jet (6 ~ 4°). Such a 
configuration is not without pre cedent and has been invoked to explain multiple breaks in 
the light curve of GRB 030329 (IBerger et al.l l2003bl ). Explaining the early emission from 
GRB080319B is beyond the scope of this work; however, we consider the implications of 
multiple-jet models in § 



3.4. GRBs 050904 and 070125 



Finally, for completeness we include here a brief summary of the primary results from 
the modeling of GRB 050904 jFrail et al.ll2~0~06h and GRB 070125 Jchandra et al.ll2008h . 



T o date, GRB 050904 is the third most distant spectroscopically confirmed GRB 



6.295; iKawai et al.ll2006l ; lHaislip et al.1 120061). The optica l and X-ray light curves exhibit a 
prominent break at t = 2.6 ± 1.0 d (jTagliaferri et al.l 120051 ) . resulting in an opening angle of 
9 pd 8°. The afterglow was notable for an extremely large inferred density: n « 700 cm -3 
for an isotropic circumburst medium. Even after applying a collimation correction, the total 
energy release from GRB 050904 was in excess of 10 52 erg, making it one of the most energetic 
explosions ever detected. 

Although Swift did not trigger immediately on GRB 070125, the 7-ray emission was 
bright enough to be detected by both Konus- Wind and RHESSI , providing superb cover- 
age of the high-energy properties of this event (IBellm et al.ll2008l ). The radio afterglow of 
GRB 070125 was one of the brightest in the Swift era, making it an ideal so urce for broad- 
band modeling. A clear break is seen in the optical light curve at t « 4d (IMirabal et al. 



20071 ). While the X-ray light curve also undergoes a steepening around this time, it occurs 
slightly later than in the optical bandpass. This may be due to the effects of inverse-Compton 
emission dominating the X-ray afterglow at this time. The circumburst density inferred for 
GRB 070125 was also relatively high (n ~ 40cm -3 ), resulting in a strongly self-absorbed 
radio spectrum. 



The relevant energy properties of all five of the events in our sample are summarized in 
Table El 
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4. Discussion 



In the previous section, we provided model fits to the broadband afterglows of five Swift 
GRBs, all of which exhibit evidence for a collimated, relativistic outflow. In some cases, the 
breaks are clearly visible across the X-ray, optical, and radio bandpasses (e.g., GRB 070125), 
while in others the data are insufficient to verify the achromatic nature of the observed break 
(e.g., the X-rays for GRB 060418). Regardless, the fact that all five events are consistent 
with a beamed outflow is in marked contrast to previous searches for jet break candidates 



in Sw ift events (jPanaitescul 1200 71 ; iKocevski &: Butlerll2008l ; iRacusin et al.ll2008bl ; iLiang et al. 
2008h . 



The most natural explanation for this discrepancy is the role of selection effects. As is 
evident from Figure [T], the heightened sensitivity of Swift is preferentially selecting GRBs 
with smaller E^^ so values compared with pie-Swift missions. Bandpass likely exacerbates this 
effect: the observed corr elation between the peak energy of the prompt 7-ray spectrum (E p ) 



and i? 7i i S o (lAmatil 120061 ) further suggests that Swift is detecting an underluminous sample 



with respect to previous missions with extended high-energy coverage. 

The result is that Swift jet breaks should occur on average later than pre- Swift events, 
making them more difficult to observe (for a given sensi tivity limit) . Such an effect was 
predicted (in the context of the structured jet model) by iPerna et all (120031 ). GRB 060418 
offers an illustrative example; with -E 7i ; S o = 10 53 erg, it falls in the 80th percentile of the Swift 



i£ 7 ,iso distribution, and with tj = 7.6 d, the break occurred at R ~ 24.5 mag. The X-rays were 
already too faint at this time to be detected, and the only reason for the optical detection 
was the deep HST imaging. Given the typical follow-up capabilities of a medium-aperture 



telescope, jet breaks are v irtually undetectable for a majority of Swift GRBs (IDai et al.ll2008 



Kocevski fc Butlerl 120081 ) . In fact, some of the faintest events (in terms of E lt is ) may be 
isotropic and still consistent with our observed energetics distribution. 



In Figure HI we plot the observed jet break times for our sample as a function of E, 



7, ISO 



Table 6. Collimation and Energetics of Swift GRBs 



GRB I tj ot E~, -Eke 





(d) 


(°) 


(10 51 erg) 


(10 51 erg) 


050820A 
050904 
060418 
070125 

080319B 


11 1+ 01 
11.1 17 

2.6 ± 1.0 
7.6+^2 
3.691*2? 
11.81?;! 


3.0 ±0.1 
8.0 

22.5±2;j> 

13.2 ± 0.6 

7 0+ - 7 
' - u -o.i 


12.9l 6 3 j 
3.0±f'f 
25.3+" 

10 - 2 ±0.i 


36. ll^ 3 

8.6±5:S 

94+ ' 22 
1 7 +0.4 
1 -' -0.2 
n or+0.38 
u - oo -0.01 
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Fig. 4. — Afterglow jet break time (rest frame) as a function of E ltiso . The dashed line 
represents a roughly constant collimation-corrected prompt energy release (there is weak 
dependence on the circumburst density and 7-ray efficiency). The events in our sample 
have on average larger jet break times than analogous pie-Swift events. This result suggests 
that the E n distribution is wider than previously th ought in a relatively mode l- independent 
manner. Pre- Swift events have been compiled from iFriedman fc Blooml (120051 ). 



compared to the pre- Swift sample from IFriedman &: Blooml (120051 ). The solid line indicates 
a constant collimation-corrected prompt energy release of -E 7i i S o ~ 10 51 erg. We wish to 



emphasize that the derived collimation angle is only weakly dependent on two model param- 
eters: n and ^ 7 (both to the 1/8 power). As a result, though the afterglow energy may be 
relatively uncertain, it is clear that the E 1 values derived for the events in our sample will 
be significantly larger than the typical pie-Swift GRB. 

The final, collimation-corrected energy release from the events in our sample, both from 
the prompt 7-ray emission and powering the afterglow, is plotted in Figure [51 Also plotted 
are analogous results from previous studies of pie-Swift afterglows (see figure caption for 
references). With the exception of the most nearby events (red diamonds), the total energy 
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release (E 1 + -Eke) from pie-Swift GRBs was relatively tightly clustered around a value of 
~ 3 x 10 51 erg (solid line in Figure [5]). Clearly the GRBs presented here are inconsistent 
with this distribution, preferentially falling at the high-energy end. Several events exceed 

10 52 erg in total e nergy release, something only achiev ed for a single event in the pre- Swift 
era (GRB 970508; hfost et al.l [20031 : berger et al]l2004l ). Much as has been seen at the low- 
energy end, our results suggest that the true energy release from GRBs is relatively broad 
and capable of extending out to at least 10 52 erg. 

Finally, we return to the issue of 7-ray efficiency As has also been seen in the pre-Swift 
era, the inferred 7-ray efficiency can often b e dramatically higher than would be predicted 
from the internal shock model (r/ 7 < 0.10; iKobayashi et al.l 1 199 71 ; iDaigne fc Mochkovitch 
19981 ). The problem may be eased to some extent by the nature of our model parameters, as 
£"ke is measured at the time of the transition from fast t o slow cooling, after which the shock 
may have lost a significant fraction of its initial energy (IChandra et al.ll2008l ). However, this 
would necessarily increase the total energy budget, possibly at times to values approaching 

10 53 erg. An alternative possibility is that th e 7-rays are produced via relativistic turbulence 
(lLazar et al.ll2009t iNarayan &: Kumarll2009l ). where the 7-ray effi ciency can approach unity 
Such a model has been already been invoked for GRB080319B (IKumar fc Narayanll2009l ). 
and clearly merits further study. 

Alternatively, if the flux distribution were dominated by small-scale fluctuations (of 
angular scale r -1 ), this may cause a large variation in the energy distribution (the "patchy 
shell" model; IKumar fc Piranll2000l ). This would have a particularly strong influence on E^, 
as the energy would be averaged over a much smaller physical scale, and could lead to large 
fluctuations in the value of r^ 7 . 



5. Conclusions 

In this work we have presented model fits to the broadband afterglows of the Swift 
GRBs050820A, 060418, and 080319B. Together with previous results from GRB 050904 and 
GRB 070125, we demonstrate all five events are consistent with our understanding of rela- 
tivistic, collimated explosions. However, the inferred opening angles for several events are 
larger than would be expected if GRBs were truly standard candles. The result is a broad 
collimation-corrected energy distribution, with some events emitting in excess of 10 52 erg. 

It is at first glance somewhat surprising that Swift has discovered a large fraction of the 
most energetic GRBs. This may be due in large part to sample size. As can be seen from 
Figure [Tj nearly all the GRBs in our sample (except GRB 060418) were brighter (in terms 
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Fig. 5. — The two dimensional (-Eke + -E-^iso) relativistic energy release from long-duration 
GRBs. Neglecting the most nearby, underluminous GRBs (GRBs 980425, 031203, and 
060218), pie-Swift events are clustered around a total energy release of p» 3 x 10 51 erg (solid 
line; dashed lines show the rms of the distribution). With the exception of GRB 060418, all 
of the events in our sample fall well outside this distribution, typically with E total > 10 52 erg. 
The solid grey line reflects a constant 7-ray efficiency of r/ 7 = 0.1. That so many of 



010222; 



Panaitescu & Kumar 




Yost et al. 


f 


2003 


): 



ble fo r the generation of the prompt 7-ray emission. Refer 
f l2002h : GRBs 990123, 99051 0, 991208, 991216 000301C, 

GRBs 970508, 98070 3 000926: iBerger et al.l J2004h : G RBs 970508, 980703: [Chevalier et al. 
fh004T): GRB 020 405: iBerg er et all (hoOlh: GRB 000418: berger et all J2003H): GRB 030329 



Soderberg et"aD d2004j): GRB 020903; iLi fc Chevalieri Jl999h : GRB 980425: ISoderberg et al. 



(l2004h : GRB 031203: ISoderberg etaU J2006bh : GRB 060218. 



of Ey teo) than nearly all of the pre-Swift GRBs having known redshifts. Even though the 
median -E 7i i SO value of Swift bursts is smaller than the pre-Swift sample, the increase in the 
total number of GRBs detected by Swift has enabled it to target the extreme ends of the 
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energy distribution. These events must be relatively rare, or they would have been easily 
detected by previous high-energy missions. 

Despite their rarity, the most energetic events provide some of the strongest const raints 
on p ossible progenitor models. The maximum energy release in magnetar models (jUsov 
19921 ) is ~ 3 x 10 52 erg, which is set by the rotational energ y of a maximally rotating neutron 



star (see, e.g., iThompson et al.l I2Q04J ; iMetzger et al.l 120071 ). The fact that several events in 



our sample approach this limit suggests that the central engine remnant should be a black 
hole (at least for these hyper-energetic GRBs). A lingering question, then, is how to produce 
the long-lived (t 3> A^crb) central engine activity seen commonly in Swift X-ray afterglows. 
While this could naturally be explained in the magnetar model, the mechanism for generating 
this activity from black holes is less clear. Though various theories have been proffered (e.g., 
Zhang et al.ll2006l ). this remains a vexing issue. 



Just as troubling is the lack of evidence for progenitor mass loss in the density profile of 
the circumburst medium of some events. In the case of radio supernovae, where the shock 
expansion is Newtonian, the circum burst medium is consistently well fit with an r~ 2 density 
profile (e.g. JSoderberg et al.ll2006al ). That this is not the case for long-du ration GRBs, given 
the preponderance of evidence for their association with massive stars (jWoosley fc Bloom 



20061 and references therein), should give pause in over-interpreting the parameters derived 
from these models. 

One possible explanation may be our relatively poor understanding of the final stages 
of stellar evolution. A variety of recent results, ranging from the discovery of fast-moving 
(v ~ 6000 km s _1 ) material ejected from rj Car (jSmithl 120081 ). to the de nse circumstellar 
material partially powering the luminous (My ~ —22 mag) SN 200 6gv ffofek et all | 2007| : 

2003; 



Smith et a 



2007), to detection of a pre-SN outburst from SN 2006jc (IPastorello et al. 



Foley et al.l 120071 ) and SN2005gl (jGal-Yam fc Leonard! 120091 ). suggest that massive stars 
undergo violent periods of episodic mass loss in the late stages of stellar evolution. Future 
theoretical progress to pin down the expected density profile surrounding massive stars in 
the latest stages of stellar evolution a t distances relevant to GRB afterglows (~ 1 pc) might 



help to resolve this discrepancy (e.g., iRamirez-Ruiz et al.ll2005l ; |Pwarkadasll2007l ). 



An additional concern, motivated by the dou ble-jet models for GRB 030329 (IBerger et al 



2003bl ) and GRB 0803 19B (IRacusin et al.ll2008al ). is our assumption that the entire relativis- 
ts outflow is collimated into a uniform jet with a single opening angle (the so-called "top-hat" 
model). A variety of other models for structured jets have been proposed, typically with 
the Lorentz f actor of the outflow varying as a function of angle from the jet axis (see, e.g., 
Granotll2007l for a review). Though the double-jet model is no doubt a bit contrived, the- 
oretical simulations of relativistic jets suggest the top-hat model is overly simplistic (e.g., 
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Zhang et al.l 12004c iTchekhovskoy et al.ll2008l ). It would also ease the sometimes extreme ef- 
ficiency requirements if the 7-ray emission were more narrowly beamed than the afterglow. 
However, such models imply the existence of so-called "on-axis" orphan afterglows, where 
the line of sight misses the 7-ray emission but an observer still sees a regular (on-axis) af- 
terglow. The lack of such events to date limits the X-ray to 7-ray beaming factor ratio to 
< 2 : 1, while future wide-field, high-cadence surveys will soon do the same for the optical 
rtNakar fc Piran]l2003h . 



Finally, we consider the future prospects for the study of the most energetic GRBs. 
As we suggested previously (IChandra et al.l 120081 ) . the recent launch of the Fermi satellite 
offers an incredible opportunity in this respect. The high-energy bandpass of the Large 
Area Telescope, extending out to hundreds of GeV, is ideally suited to target the bright 
end of the -E 7i i SO distribution. Already in less th an a year of operation, Fermi has det ected 



two of the br ightest GRBs eve r GRB080916C (Abdq et all 12009: iGreiner et al J 1200a ) and 



GRB 090323 flOhno et ajJIgOOJ ; iGolenetskii et al. 



20091 ). with -E 7iiso in excess of 5 x 10 54 erg. 



Coupled with the large Lorentz factor required to produce GeV photons, multi-wavelength 
campaigns targeted at such events are well positioned to search for early jet breaks when the 
afterglow is still bright. Together, synergistic Swift and Fermi observations in the coming 
years should be able to shed a good deal of light on the high end of the GRB energy 
distribution. 
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• 8.5 GHz 

— ISM-1 : 9=6.6° 
---ISM-1: Isotropic 




Time Since Burst (d) 



Fig. 6. — The X-ray (previous page, top panel), optical (previous page, bottom panel), and 
radio (this page) light curves of GRB 050820A. Radio observations at t < 5 d a re left out of 



the mo del, as the emission at this time is likely dominated by the reverse shock (ICenko et al. 



2006bJ). The preferred model (ISM-1; Table [3]) is plotted as a solid line, while the identical 



model for an isotropic explosion is shown as a dashed line. The optical data have been scaled 
by the model flux to match the i?c-band. Both the X-ray and optical bandpasses show a 
clear break at t ~ 10 d. The radio is not sufficient to distinguish between an isotropic and a 
collimated outflow. 
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Fig. 7. — The X-ray (previous page, top panel), optical (previous page, bottom panel), and 
radio (this page) light curves of GRB 060418. The preferred model (Wind-Table HJ) is plotted 
as a solid line, while the identical model for an isotropic explosion is shown as a dashed line. 
The optical data have been scaled by the model flux to match the i?c-band. The optical 
shows a clear break at t ~ 10 d, which is also favored in the radio. The X-ray afterglow is 
too faint to be detected at this time. 
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• X-ray(2keV) 
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Time Since Burst (d) 



Fig. 8. — The X-ray (previous page, top panel), optical (previous page, bottom panel), and 
radio (this page) light curves of GRB080319B. The preferred model (Wind-2-Table H|) is 
plotted as a solid line, while the identical model for an isotropic explosion is shown as a 
dashed line. Both the X-ray and optical show a break at t ~ 10 d. The radio afterglow is too 
faint at this time to see evidence for collimation. We have left out the late-time (t > 10 d) 
i', and d data from our fits due to the presence of contaminating SN emission. 



Table 7. Optical/NIR Observations of GRB 060418 



UT Date a Telescope / Instrument Time Since Burst b Filter Exposure Time Magnitude 
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Rc 


1800.0 


> 21.02 
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Table 7 — Continued 



UT Date a 


Telescope / Instrument 


Time Since Burst b 
00 


Filter 


Exposure Time 
(s) 


Magnitude 


2006 Apr 19.3581 


PAIRITEL 


1.079 


X 


10 s 


J 


2283.8 


> 19.06 


2006 Apr 19.3581 


PAIRITEL 


1.079 


X 


10 s 


H 


2283.8 


> 18.61 


2006 Apr 19.3581 


PAIRITEL 


1.079 


X 


10 s 


K s 


2283.8 


> 17.89 


2006 Apr 19.3669 


P60 


1.081 


X 


10 6 


i' 


1800.0 


> 21.49 


2006 Apr 19.3782 


P60 


1.091 


X 


10 5 


V C 


1800.0 


> 21.36 


2006 Apr 20.3405 


P60 


1.920 


X 


10 s 


Rc 


1800.0 


> 21.75 


2006 Apr 20.3535 


P60 


1.931 


X 


10 s 


i! 


1800.0 


> 22.19 


2006 Apr 19.2870 


P200 / LFC 


1.002 


X 


10 5 




300.0 


22.07 ±0.10 


2006 Apr 19.2954 


P200 / LFC 


1.009 


X 


10 s 


i' 


300.0 


21.73 ±0.10 


2006 Apr 19.3033 


P200 / LFC 


1.016 


X 


10 s 


t! 


300.0 


21.52 ±0.29 


2006 Apr 19.3122 


P200 / LFC 


1.024 


X 


10 s 


</ 


300.0 


22.75 ±0.15 


2006 Apr 22.5584 


Keck / LRIS 


3.833 


X 


10 s 


Rc 


1200.0 


23.31 ±0.08 


2006 Apr 22.5590 


Keck / LRIS 


3.835 


X 


10 s 


</ 


1500.0 


> 23.08 


2006 May 9.4161 


HST 1 ACS 


1.839 


X 


10 e 


F625W 


4220.0 


27.06 ±0.15 


2006 May 9.6130 


HST 1 ACS 


1.856 


X 


10 e 


F775W 


4220.0 


26.96 ±0.19 


2006 May 20.2953 


HST 1 ACS 


2.779 


X 


10" 


F625W 


5500.0 


> 27.53 


2006 May 20.4926 


HST 1 ACS 


2.796 


X 


10" 


F775W 


3700.0 


> 26.93 


2006 Jun 1.7528 


HST j ACS 


3.855 


X 


10 6 


F625W 


8772.0 


> 28.04 


2006 Jun 2.0193 


HST 1 ACS 


3.879 


X 


10 e 


F775W 


8772.0 


> 27.51 


2006 Jul 11.1211 


HST 1 ACS 


7.343 


X 


10 6 


F625W 


8772.0 


> 28.16 


2006 Jul 12.5207 


HST 1 ACS 


7.378 


X 


10 6 


F775W 


8772.0 


> 28.10 


2006 Jul 13.1036 


HST 1 ACS 


7.428 


X 


10 6 


F555W 


4386.0 


> 27.53 



a UT at beginning of exposure. 

b Time from midpoint of exposure to Swift-BAT trigger. 

c R eported magnitudes have not been corrected for Galactic extinction (E(B — V) = 0.22 mag; Schlegel etalj 
Il998h . Observation s in the Vc,Rc, J t H, and K s niters are referenced to Vega, while all other niters use the AB 
magnitude system <Oke fc Gunnlll983t) . 



